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1 Introduction 

Let f2 be a bounded domain of R"(n = 1,2,3) with smooth boundary. Let T be an 
arbitrary positive number. Let u = u{x,t) satisfy 

ut = Au + f{x,t) mQx]0, T[. 

In this paper we consider an inverse problem for the heat equation. The problem is 
Inverse Problem. Assume that there exist a non negative number Tq less than T and 
point Xq E fl such that /(xq, Tq) ^ and f{x, t) = for all < t < Tq and all x E fl. 
Extract Tq and information about the set {x E Q \ f^x.To) ^ 0} from the data M|anx]o,r[5 
9M/(9i^|anx]o,T[ and u{ ■ ,0). 

The number Tq and the set {x E VL \ /(x, Tq) 7^ 0} are the time and position when and 
where the heat source /(x, t) firstly appeared. Note that one may assume that ■ , 0) = 
in VL and in what follows we do so. 
Some related known results should be mentioned. 

Yamatani-Ohnaka [23] considered the case when the source f{x,t) takes the form 



N 

fi.x,t) = J2pj^(^ ~ Xj^t- tj, 
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where the points Xj e times tj G ]0, T[, strength Pj{< 0) and N are all unknown. 
Note that, in this case Tq = mmj{tj}, formally the set {x E ^ Tq) 7^ 0} coincides 
with {xj I To = tj} and the source becomes inactive after the time max^ tj . They made 
use of special solutions of the backward heat equation Vt + Av = in fix ]0, T'[ where 
T >T' > maxj tj having a large parameter c. Those solutions have singularity at t = T' 
on the given plane outside Q and are vanishing at i t T" on Q. Those are constructed from 
the fundamental solution of the backward heat equation in one-dimensional space. Then, 
using integration by parts, they obtain a system of equations involving those solutions, 
unknown Xj,tj,pj and N. Carefully analyzing the asymptotic behaviour of the system as 
c — > 00, they gave a reconstruction formula of the source itself. 

Yamamoto [22] considered an inverse source problem for the wave equation, however, 
it is possible to apply his method to the inverse source problem for the heat equation. 
One interpretation of his method for the heat equation is the following. The assumption 
on the heat source is that 

f{x,t) = ait)fix) 

where a is a known function with o"(0) 7^ and / is unknown. From our point of view 
this is the case when Tq is known and Tq = 0. 
Then, his method consists of two parts: 

• a way of calculating the data dw/du\dnx ]0, T[ where w solves 

wt = Aw infix ]0, T[, 
w{x, 0) = f{x) infi, 

w{x,t) = on^fixjO, T[, 

from the data wjanx ]0, T[, du/di'\dQx]0, T[ and the time derivative {d / dt)du/ dv\dQ,x]Q, Tf; 

• reconstruction formula of f{x) itself from the calculated data dw/dulg^x ]0, T[. 
The first part can be done by solving a Volterra equation of the second kind. For es- 
tablishing the second part the null- controllability for the heat equation [7, 20] and the 
completeness of the eigenfunctions of the Dirichlet Laplacian in fi are essential. 

In [3] El Badia-Ha Duong also considered the case when Tq = 0. They assume that: 
the source is supported in a subset of fi (mainly the point sources) that is independent of 
the time, the strength of the source on the set does not depend on the position and has 
a definite signature in an interval ]0, with < T; the source becomes inactive after 
the time T^. 

Their method also consists of two parts: 

• a way of calculating the temperature in fi at t = Ti from the data uIq^x ]Ti, T[ and 
du/duldfiX ]Ti, T[ where Ti is an arbitrary fixed time with T^, < Ti < T; 

• extracting moments of the source from u{x,Ti) and the data and determination of 
the sources form the moments. 

The first part is also an application of the null-controllability for the heat equation and the 
completeness of the eigenfunctions of the Dirichlet Laplacian in fi. In the second part they 
start with the determination of the number and location of point sources. This is purely 
algebraic and can be done by using a previous result in the inverse point source problem 
for eUiptic equation [2] . The next step is the determination of the strength of the point 
sources from some exponential moments which is a combination of the Miintz's theorem 
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(e.g., [8]) and an invertibility of a matrix involving special solutions of the Helmholtz 
equations. As they pointed out the choice of the matrix is not constructive. These steps 
can be done by a combination of an integration by parts and special solutions for the 
backward heat equation. 

It should be pointed out that the controllability for the heat equation plays an im- 
portant role also in the problem of recovering the conductivity coefficient of the heat 
equation, see Avdonin-Behshev-Rozkov [1]. 

In this paper we apply the enclosure method to Inverse Problem which was intro- 
duced by Ikehata [13] for inverse boundary value problems for the elliptic equations. In [10] 
Ikehata considered in two-dimensions, an inverse source problem for the inhomogeneous 
Helmholtz equation Au + k^u — f with the unknown source / having the form xd{x)p{x) 
and established an extraction formula of the convex hull of D provided both D and p 
are unknown and D is polygonal. Now the enclosure method has been applied to inverse 
boundary value problems [11, 12, 18], a Cauchy problem for the stationary Schrodinger 
equation [14] and inverse obstacle scattering problems [15, 16, 17]. The common character 
of these problems is: the governing equation is elliptic. 

It is quite curious to consider whether one can apply the idea of the enclosure method 
to the case when the governing equation is non elliptic since there are many inverse 
problems whose governing equations are non elliptic. The inverse source problem for the 
heat equation is a typical and important one. The aim of this paper is to introduce a 
direct approach to the inverse source problem for the heat equation by employing the idea 
of the enclosure method. 
The main feature of our method is: 

• it is based on a simple one line formula; 

• we do not make use of the exact controllability of the heat equation nor the com- 
pleteness of the eigenfunctions of the Dirichlet Laplacian in Q; 

• the assumptions on the unknown source is quite general; 

• the method provides us a brief information about the time and the location when 
and where the unknown source firstly appeared instead of the detailed information of the 
source. 

For several other formulations of the inverse source problem for the heat equation and 
results including uniqueness and stability see Cannon- Esteva [5] , Cannon-DuChateau [4] , 
Isakov [19], HetUich-Rundell [9], recent pubhcation of Trong- Long- Alain [21] and refer- 
ences therein. 

2 Statement of results 

The results of this paper are divided into two parts. The first one is a formula for 
extracting Tq only and second is for extracting both Tq and information about the location 
of the set {x E Q \f{x, Tq) ^ 0}. 
2.1. Extracting Tq. 

The list of the assumptions on the source f{x,t) is the following. 

(A.l) f{x,t) takes the form 

f{x,t) = XD{x,t)p{x,t) 
where D C Q x [Tq, T] is a Lebesgue measurable set; p is essentially bounded on D. 
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Thus the source may appear at the time Tq firstly and the points in D n (R" x {Tq}). 
The next assumption describes the standing behaviour of the source. 

(A. 2) There exist positive numbers 5 < T — Tq, Ci, C2 and p G [0, oo[ such that the 
n-dimensional Lebesgue measure of the set D[s) = D (R" x {Tq + s}) has the lower 
estimate: 

\D{s)\> CisP for almost alls e [0, 61]; 
the strength of the source p satisfies 

p{x, t) > C2 for almost all {x, t) e D n x [Tq, Tq + 5]) 

or 

-p{x, t) > C2 for almost all {x, t) e Dn{W x [Tq, Tq + 5]). 

Definition 2.1. Given u; e -S'""-'^, s e R define the indicator function I^{t;s) by the 
formula 

Ur; s) = J^^ J^^ (^^u - dSdt, r > 

where v = v{x,t) = e^^'^~'^^. 

The function v takes positive values and satisfies the backward heat equation Vt + Av — 
in the whole space. Moreover e'^^'w has the special character: 

• if t > s, then lim^ >oo e'^^v{x, t) = 0; 

• if t < s, then lim^ e^^v{x, t) — 00. 

Note that the function w{x) — e^^''^ satisfies the equation 

Aw — Tw = 

in the whole space. Since r > 0, this is not the Helmholtz equation. In the following we 

make use of the asymptotic behaviour of w as r — > 00. 

Theorem 2.1. Assume that assumptions (A.l), (A. 2) are satisfied. As r — > 00 we have 

hm i^st^IlM ^ (2.1) 

r — >oo f 

Moreover, we have: 

if s < Tq, then lim^ = 0; 

if s > Tq, then lim^ »oo \Iui{t]si)\ — 00. 

The result can be generalized to the case when the governing equation has variable 
coefficients. In that case there is an interesting relationship between the high frequency 
asymptotic solution of the corresponding wave equation [6] and that of the formal adjoint 
equation of the heat equation. This is clarified in Section 3. 

2.2. Extracting both Tq and information about the set {x e | /(x, Tq) 7^ 0}. 

This subsection starts with the following observation. 
Let n > 1. Let c 7^ be an arbitrary number. Let 

= e-(^-^V-^ (2.2) 
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where r satisfies r > c ^ and 



CT \ -\- —ijj I , n > 2, 



CT I l + 1 - ^ I , n= 1. 



The function f is a complex valued function and satisfies the backward heat equation 
Vt + Av = 0. Moreover e'^'^v has the special character: 

• if s < t — cx ■ uj, then lim^ ,00 e.'^^\v{x, t)\ =0; 

• ifs>t — cx-cu, then liniT- >oo e'^'^\v{x, t)\ = 00. 

So one can expect that: using this function, one may obtain more information about the 
location and shape of D in the space time. We show that the method [14] is applicable 
in the case when n = 1,2 (at least). 

Definition 2.2. Let n > 2. Given c > 0, s G R, a;, a;^ G 5""^ with a; ■ a;^ = define 
the indicator function Iw,w^,c{^') s) by the formula 

I.,.^Ar; s) = e- I i^-u - -vj dSdt, r > 
where v is the function given by (2.2). 

Let n — 1. Given c 0, s e R, define the indicator function Ic{t; s) by the formula 

where v is the function given by (2.2). 

Let n>2. Given u; e S''-'^ and c> define the unit vector u;{c) in R"+^ = R" x R 
directed into the half space i < by the formula 

Then 

t — cx ■ uj — —V + 1 ^ ^ ■ u>{c). 

Thus we see that: 

• if — s/ a/c^ + 1 > (x tY ■ ^(c), then lim^- >oo e'^^\v{x, t)\ =0; 

• if — s/vc^ + l < [x t)'^ ■ uj{c), then then lim,- ^^e'^^\v{x,t)\ — 00; 

First we consider the case when n = 2. 

We assume that the unknown source takes the form 

N 

fi^^t) = Y,XPjx[Tj,T]{x,t)pj{x,t) 

where 

• Pj CQ is given by the interior of a polygon and Tj satisfies < 7} < T; 
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• iij ^ f, then Pj n Pf = 0; 

• for each j pj e L°°{PjX]Tj, T[) and for each vertices p of the convex huU of Pj pj 
coincides with a Holder continuous function with exponent 9j e]0, 1] in a neighbourhood 
of {p,Tj) that does not vanish at {p,Tj). 

Wc set D = uf^-^{PjX]Tj, T[). 

Theorem 2.2. Let a;(c) be regular with respect to D. Assume that 

sup ( ^Yuj{c)<hD{uj{c)). (2.3) 

Then the formula 

y log 14,0;^ ,c(t; 0)1 r^—-rr, / i 

hm ■ — '■ = Vc^ + 1 /iD c , 

is valid. 

Moreover, we have: 



ifkniuj^c)) < -si a /c^ + I , then hm^^oo \Iui,u,^,cij; s)\ = 0; 
ifhoioJ^c)) > -s/Vc^ + 1; then HuIt—^oo \Ilo,uj^,c{'^'i ^)\ = 

The condition (2.3) means that the set Q x {T} is located in the half space (x, t)^ •a;(c) < 
hD{i^{c)). This can be satisfied in the case when: T is given and c is sufficiently small; c 
is given and T is sufficiently large. 

Next we consider the case when n = 1. 
We assume that the unknown source takes the form 

N 

fi^^t) = J2XPjX[Tj,T]{x,t)pj{x,t) 
j=l 

where 

• Pj C is given by the interior of an interval and Tj satisfies < 7} < T; 

• if f, then Pj n Pj/ = 0; 

• for each j pj G L°°{PjX]Tj,T[) and for each end point p of Pj pj coincides with a 
function of class in a neighbourhood of the point (p, Tj) that does not vanish at (p, Tj). 
We set D = uf^^{PjX]Tj, T[). 

Given c 7^ define the unit vector u;{c) in R^+^ = R x R directed into the half space 
i < by the formula 

Theorem 2.3. Let uj{c) he regular with respect to D. Assume that 

sup ( ^ ) ■uj{c) < hD{uj{c)). 



Then the formula 



is valid. 

Moreover, we have: 



log|/c(r;0)| Ay— T, / / 
hm — V c^ + 1 hoii^ic)), 
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if hD{^'{cj) < —sj^/c? + 1, then lim^ |-^c(t; s)\ = 0; 

if hD{uj{c)) > —s/y/c^ + 1, then IrniT >oo |-^c('?';s)| = oo. 

If Q C R^, the situation seems more complicated. We will consider the case in the future 
work. 

3 Proof of Theorem 2.1 and its generalization 

3.1. Proof of Theorem 2.1. 

Since we have the trivial identity 

7,(T;s) = e-(^-^°)/,(T;ro) Vs, (3.1) 

it suffices to study the asymptotic behaviour of the indicator function at s = Tq. Integra- 
tion by parts gives 

/ / f(x,t)v(x,t)dxdt— / u(x,T)v(x,T)dx + / / (—u — —v)dSdt. 
Jo Jn Jn Jo Jdct ou ou 

Thus we have the representation of the indicator function at s = Tq: 

I^r; To) = e^^° r [ f{x, t)v{x, t)dxdt - e^^° / u{x, T)v{x, T)dx. (3.2) 
JO Jn Jn 

Since T > Tq, it is easy to see that the absolute value of the second term of the right 
hand side is dominated by 

/ k(x,r)|dxe^^"P-6n^--'e-^(^-^°) = 0(e-^(^-^°)/2) (3.3) 
Jn 

as r — > oo. 

On the other hand, from the assumptions (A.l) and (A. 2) one knows that the first term 
of the right hand side of (3.2) takes the form 



e^^° / p(x,t)v(x,t)dxdt 
Jd 



gTTo 



/ p(x,t)e^''-^-^^dxdt + e^^° [ p(x,t)e^''-^-^^dxdt. 

Jon {R"x [To,To+S]) Jd\(R"x [To,To+S]) 

(3.4) 

Since the set D \ (R" x [Tq, Tq + S] is contained in R" x [Tq + S, T], we have 

e^^» / p(x,t)e^''-''-^*dxdt ^ 0(e-^^/'^) (3.5) 

iD\(R-"x[ro,ro+<5] 

as T — > oo. 

Here we prepare 

Lemma 3.1. As r — > oo we have 

irie^^V-(f+') < e^^°| / pix,t)e^^-^-^'dxdt\ < K^e^^' (3.6) 

Jon (R"x [To,To+<5]) 

where Ki, K2, K^, K4 are constant and K^^K^ > 0. 
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Proof. 

We just describe the key point of the proof of the left half of the inequahties in the 
case when 

p{x, t) > C2 for almost all {x, t) e D n {W x [Tq, Tq + S]). 

We have 



/ pix,t)e^^-^-^'dxdt>C2e'''^° j e^^-^-^'dxdt 

JDniR^x [To,To+S]) ■/I3n(R"x [To,To+<5]) 

= C2 r°^^ e-^(*-^°) I / e^^-^dx 1 dt 

JTo \JD{t-To) J 

= C2 t e-^' { I ev^^-^dx I ds 

Jo \JD{s) J 

Jt inf a; ■ a; r& 
>C2e x^^ / \D{s)\e-^'ds 
Jo 

\/t inf a; ■ a; 
>CiC2e / s^e-^'ds 

Jo 

Jo 

□ 

Now Theorem 2.1 is a consequence of (3.1) to (3.6). 
3.2. A generalization 

In this subsection we give a generalization of Theorem 2.1 to the heat equation with 
variable coefficients: 

aut^V -iSIu + Fix.t) 

where a = a{x) is a smooth function with positive values and 7 = 7(2;) is a n x n-real 
symmetric positive definite matrix valued smooth function. We apply the idea of the 
enclosure method to this equation. For the purpose we construct a solution with large 
parameter r of the equation 

avt + V ■ ^Vv = (3.7) 

which plays the role of the function e'^^'^"'^* in Theorem 2.1. 

First we consider how to construct a solution of the equation (3.7) in the form: 

v{x, t) = e~'^^w{x). 

Substituting this into (3.7), we have 

V • 7VW - Taw = 0. (3.8) 

Note that r > unlike the case when the governing equation is the wave equation ([6]). 
Let be a smooth function. The change of the dependent variable 

w — e^'^w 
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gives 

V • jWw' + (27V(^ • Vw' + (V • jWip)w') + T{{jWip ■ Wip) -a)w' = 0. 
Thus if (f satisfies the eikonal equation 

7V(^ ■V(p = a, (3.9) 

then w' satisfies 

V • 7Vw' + ^/t {2'yVip ■ Vw' + (V • ^Vip)w') = 0. (3.10) 

The next lemma is crucial for the construction of the exact solution from the asymptotic 
solution. 

Proposition 3.2. Let r be an arbitrary positive number. Given f £ L^(Q) there exists a 
unique weak solution u e H^{VL) of the elliptic problem 

V • -fVu + ^/t (27V(/7 • Vit + (V • -iV^p)u) ^ f inn, 
u — ondfl. 

Moreover u has the estimate 

< C||/|U2(n) (3.11) 
where C is a positive constant independent of t. 

Proof. Using a change of dependent variable, one can easily deduce the uniqueness and 
existence of the solution for the corresponding fact for the equation (3.8) in Q,. The 
positivity of r is essential. The problem is the estimate (3.11). Since 

(27V¥? ■ Vu)u = -fVif ■ V(m^), 

integration by parts gives 

(27 • + (V • 'yVip)u) udx = 0. 
Then from the equation for u we obtain 

/ 7Vti • Vudx = — fudx. 

JQ JQ 

A standard argument gives (3.11). 
□ 

From the equation one knows that the solution u in Proposition 3.2 satisfies V ■ ^Vu G 
L^(f2). A combination of the standard elliptic estimate for the operator V • 7V and (3.11) 
yields 

\\u\\HHn)<C'{^+l)\\f\\L2^a). 
Then the Sobolev imbedding yields 

IkllLoc(n) < C"{V^+l)\\f\\mn). (3.12) 
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Given > 1 set 

where ao, • • • , oat satisfy the transport equations 

2^V(p ■ Vao + (V • 7V<^)ao = 0, 
■ Vaj+i + (V • 7V<^)%+i = -V • -fVaj, j = 0, • • • , AT - i. 
Then we see that 

V • 7Vaiv(x) 



(3.13) 



AT 



This means that is an asymptotic solution of the equation (3.10). Then Proposition 
3.2 ensures the unique solvabihty of the problem 

V • 7Vm + Vt (27V<^ • Vm + (V • 7V(^)m) = ' /^[^^ in 

u = ondfl. 

Write u — Rn- Then the function 

is an exact solution of the equation (2.10) and, by virtue of the estimate (3.11), we have 

\\RN\\H^{n) = O {j-;j^y^ ■ 



Thus we have 

w 



Note that w"^ is the unique solution of the problem 

V • 7V'u + ^/t (27V(^ • Vm + (V • 7V(^)'u) = in O, 

u — w'f^ ondfl. 

Summing up, we have 

Theorem 3.3. Let if satisfy the eikonal equation (3.9). Given N > 1 let the functions 
ao, - ■ ■ ,aN satisfy the transport equations (3.13). Then the unique solution of the elliptic 
problem 

V • 7Vu' — Taw = mJl, 

^ a, 



^ = e^^V-^ ondn 
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satisfies, as r — oo 



N-l 



iv^y 



O 



N 



Note that, from (3.12) we also have 



where > 2. 

Thus we have a special solution v of the backward heat equation (3.7) in Q x R+ 
having the form 

„(M) = e-^'e-{|'gM.0(^)) (3.15) 

as r — > oo provided Lp, a^, ■ ■ ■ , a^v are given. 

Now we define another indicator function. 
Definition 3.1. Given s e R define the indicator function Itpir; s) by the formula 

where v is given by (3.15). 

For the description of the result we introduce the assumption instead of (A. 2): 

(A. 2)' There exist positive numbers 5 < T — Tq, Ci, C2 and p G [0, oo[ such that the 
n-dimensional Lebesgue measure of the set D{s) = D (R" x {Tq + s}) has the lower 
estimate: 

\D{s)\ > Cis^ for almost alls e [0, Si\; 
the strength of the source p multiplied by ao{x) satisfies 

p{x, t)ao{x) > C2 for almost all {x, t) e D n (R" x [Tq, Tq + 6]) 

or 

-p{x, t)ao{x) > C2 for almost aU {x, t) e D n (R" x [Tq, Tq + d]). 

We obtain 

Theorem 3.4. Assume that the assumptions (A.l), (A. 2)' are satisfied. Then, as r — 
00 we have 



T ►OO J- 

Moreover, we have: 

if s < Tq, then lim.^ >oo \I<p{t] s) \ = 0; 

if s > Tq, then lim^^ ►oo \I<fi{'T] s)\ = 00. 

Proof. (3.12) for A'" 2 gives the estimate 



|w -ao(-)IU°°(f2) = 
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Then, it follows from (A. 2)' that there exists tq > such that, for all r > tq 
p{x, t)w{x) > C2/2 for almost all {x, t) e D n (R" x [Tq, Tq + 6]) 

or 

-p{x, t)w{x) > C2/2 for almost all {x, t) e D D {W x [To, Tq + S]). 

Hereafter we take the same course as the proof of Theorem 2.1. 

□ 

Example 3.1. Consider the simplest case: a = 1 and 'j = In- 
Given an arbitrary point p G R" \ Cl the function Lp given by the formula 

(p{x;p) — ±\x — p\, X E R" \ {p} 

satisfies the equation (3.9) for a and 7 specified above. Since 

, n — 1 

A(/P = ±l r, 

\x — p\ 

the transport equation for oq in (3.12) becomes 

Tl — 1 

{x -p) ■ Voo H = 

and the method of bicharacteristics yields 

ao{x) = ao(^^ +p)\x- pj-^^-^V^ 
\x — p\ 

Thus if we specify the boundary value of Qq on the unit sphere centered at p as 

QqIx) = 1 on |a; — p| = 1, 

then we have 

One can also obtain ai explicitly. The point is: the function oq specified above never 
vanishes on Q. From Theorem 3.3 one obtains the solution v — v±{x, t;p) of the backward 
heat equation having the form 

as T — > 00. 

4 Proofs of Theorems 2.2 and 2.3 

4.1. Asymptotic behaviour of an integral 

In this subsection, we consider the case when n > 2. we assume that D C R"+^ is a 
finite cone with a vertex at p and a bottom face Q. More precisely, D takes the form 

g 

D = Uo<s<s{p +-^{y-p)\y&Q} 
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and Q is a bounded open subset of the n-dimensional hyper plane 

j ■ ^(c) = P ■ uj{c) - 5 

where 5 is a positive number. 

We consider the integral of the function v{x,t) given by (2.2) over D: 



D 



I{t) = l_v{x,t)dxdt 

Lemma 4.1. The limit 



lim — (cr)"'*'^ exp |— Vc^ + irp • a;(c)| exp |— icry^ 1 — -^P ■ ^ ) |^^^'' 



exists and has the integral representation: 



= 25 / , ^ ''^'^ . ,,n.^ ■ (4.1) 

-liy-p) ■ 



c 







Proof. 
Write 



exp j-Vc^ + Irp • a;(c)| exp |-icr^l - -^p ' (^^q j |-^(''') 

= /^exp ( ( n - ^^'^ ^^p^^^i - i ( ( t ) " ( ""o ) 



dxdt 







Then, a combination of Lebesgue's dominated convergence theorem and the formula 

/ e"e"«e'"«rfe = T. -wT' « ^ R 

Jo (l-za)"+i 
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gives 



n! r dS{y) 

2 

□ 

4.2. Non vanishing of 

It seems difficult to check Kd 7^ by using (4.1) directly. In this subsection, we give 
an alternative expression of Kd which yields Kd 7^ directly in the case when n — 2 

provided D C R^^^ is given by the interior of a tetrahedron with a vertex at P with the 
bottom face Q located on the plane (x, ty ■ uj{c) = p ■ uj{c) — 6. The set Q is given by the 
interior of the triangle with vertices yi, 1/2, 1/3 which satisfy yj ■ a;(c) = p ■ a;(c) — 6. 
Set 

A = {{a, P)\a,P>0,a + P<l}. 
Define the maps (pi,(p2:V3 '■ ^ — ^ R^"*"^ by the formulae 

(pi{a, (3) =p + a{yi - p) + (5{y2 - p), 
ip2{a, (3)=p + a{y2 - p) + (3{yz - p), 
ips{a, (3)^p + a(ys - p) + (3{yi - p) 

and set 

A,- = ¥^,-(A). 

Renumbering yi,y2,y3 if necessary, one has the decomposition of dD: 

ai:> = Q u Ai u A2 u A3. 

Q and A^-, j = 1, 2, 3 satisfy Q U Aj = and Int Aj n Int Aj, = if j 7^ /. 

Let I' be the unit outward normal vector field to dD. Since u takes a constant vector on 

each Aj, we denote the vector by Uj. For simplicity of description we identify 3 + 1 with 
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1 and 1 — 1 with 3. Renumbering yi,y2, Vs if necessary, one may assume, in advance, that 

(i/i X 1/3) •a;(c) < 0, 

(u2 X ui) ■ u{c) < 0, 

(i/3 X i/2) • ou{c) < 0. 

Proposition 4.2. T/ie formula 
is valid where 

Proof. Let a be an arbitrary constant complex vector in four-dimensions. Since 

y{x,t) ■ {va) = I I J ■ a U, 



we have 



( ^ I vdxdt — ^a ■ Vj / vdS{y) — a ■ uj{c) / vdS{y). (4.3) 

It is easy to see that, as r — 00 

/ = O (e^^^(f'^(^)-^)) . (4.4) 

On the other hand, using a similar computation as in [14] , we have 

vdS{y) 

(jj^ \ 1 5^\{iyj X X (//j+i X i^j] 



exp |V + Irp • a;(c)| exp < icrW 1 • 



X 



y exp | — \/c^ + lT(a + /5)(5|exp |— icry^ 1 ^(5(q;ci + /5c2) |(iQ;(i/? 



(4.5) 
where 

{uj X • ( '^Q J i^j+i X t^i) • ( '^0 ) 



Cl — p s . . , C2 



{uj X Uj^i) ■ u;{c) ' (i/j+i X Uj) ■ u;{c) 
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A combination of the change of variables and Lebesgue's dominated convergence the- 
orem gives 

(5t)^ J exp |— Vc^ + 1t(q; + P)S^ exp | —icT^ 1 —5{aci + (3c2) \dad(5 



- J df3 J da exp |— Vc^ + l{a + /3) | exp | — ^cy 1 ^ (aci + /3c2) 

— )• y d(3 J dctexp |— V + 1(q; + exp {—ic{aci + /3c2)} 

/•cx) /-oo f C 1 

= (c^ + l)"^y^ daexp{-(Q; + /3)}exp|-i--^===(Q;ci + /3c2)| ^4_g^ 



= (c^ + l)-^- 



1 

(Vc2 + 1 + icci) (Vc2 + 1 + icC2) 

^ i^j X ^j-i) • ^'(c) {uj+i X z/j) ■ a;(c) 

A combination of (4.5) and (4.6) gives 

exp j-Vc^ + Irp • a;(c) I exp | -icT^ 1 - — p ' ^ j | ^ . ""^^^y^ 



Then from this, (4.3) and (4.4) we have 



X exp j-Vc^ + Irp ■ a;(c)| exp j-icr^ 1 - -^p ' j ^-^W (4-7) 



E 



t (i^j X i/j_i) • (l/^+i X Uj) ■ ■§ 
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A combination of Lemma 4.1 and (4.7) yields 



= lim <- + Vl-^-") -a 



2 

X 

2 



i(cr)^exp |-Vc2 + lTp-a;(c)|exp |-icy^l - ^^-p ■ ^ ) j"^'"^'' 



3- 



1x (i^i X i/j_l) •i?(i/j+l X i/j) -i? 



Since a is arbitrary, we obtain (4.2). 

□ 

Since V\.V2. v-?, arc linearly independent, one has the conclusion 
Corollary 4.3. ^ jor all c> 0. 

4.3. Proof of Theorem 2.2 

Integration by parts gives 

Lu;^c{t;0)= r I f(x,t)v(x,t)dxdt- [ u(x,T)v(x,T)dx. 
Jo Jn Jci 

Then (4.8) yields that the integral 

^-rV^hn{u,{c)) f u(x,T)v(x,T)dx 

Jn 



is exponentially decaying as r — > oo. Thus the integral e '^^'^^'^^^'^^'^^'^^^ Iaj,u}-^,c(T", 0) 
modulo exponentially decaying as r — > oo, coincides with the integral 

^-rV^hniu^ic)) r f f(^r,^t)v{x,t)dx 

Jo Jn 

Since uj{c) is regular with respect to D, there exits a unique point p on dD such that 
p ■ cj(c) = /;,£)(u;(c)). This p should belong to some x {Tj-g} since the time component 
of uj{c) is negative. Then PjiX ]Tj/, T[ with f ^ Jq should be contained in the half space 
{x, ty ■uj{c) < hoi^jj) — S with a positive number 5 independent of /. Thus the right hand 
side of (4.9) modulo exponentially decaying as r — > oo, coincides with 



hoi^ic)) [p.^[x,t)e-^'e'''dxdt (4.10) 

J D 



where 

D' = {Pj,x]Tj, T[) n {{x,t) I {x,tf ■u>hD{oj)- 5}. 
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Choosing a smaller 5 if necessary, one may assume that this D' is given by the interior of 
a tetrahedron with vertex p and has the factorization as that of D in subsection 4.1: 

D' = Uo<s<s{p +^{y-p)\y ^Q} 

where Q is the interior of a triangle lying on the plane {x, t)'^ ■ ui{c) = /iD(a;) — 5. This is 
the key point of the proof. Divide (4.10) into two parts: 

Jd' Jd' 
= 1 + 11. 

From a combination of Lemma 4.1 and Corollary 4.3 one knows that / as r — > oo decays 
really algebraically: 

\im Jcrf\I\ = \Kd'\\PjM\- 



T >00 



Using the Holder continuity of pj^, one can also show easily 

T'\ii\^o{T-'^oy 

Summing up, one concludes that, as — > oo 

Then all of the statements follow from this asymptotic formula. This completes the proof 

of Theorem 2.2. 

4.4. Proof of Theorem 2.3 

This is the case when n — 1. Since 



z ■ z = T + i2cV^Jl 



9 ' 



we have 



v(x, t) = exp jVc^ + It ^ ^ j 



^(c) y exp ^cr ( M . ( J Jl - 



Note that this function is oscillatory higher than the case when n > 2 because of the 
existence of the growing factor in the imaginary part of the phase function. 
As same as the proof of Theorem 2.2 it suffices to show that the integral 

exp |— V + 1'7'/id(<^(c))| J v{x,t)p{x,t)dxdt 

is really algebraic decaying as r — > oo. 
Remarks are in order. 

• One may assume that D is given by the interior of a triangle with vertices p, yo and 
yi with p ■ uj — hD{oj{c)), yo ■ a;(c) = yi ■ a;(c) = /id(<^(c)) — 5 and has the form 

g 

D = Uo<s<s{p +-^{y-p)\y&Q} 
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where Q is the segment with endpoints yo, yi. And also one may assume that pj,j ~ 
,N are given by a single p = pj^ for some jo. 

• One may assume that the vectors yi—p and yo —p coincide with \yi — p|(0 1)-^ and 
\yo — sgnc 0)"^, respectively. 

• Since we have assumed that p coincides with a function of class in a neighbour- 
hood of the point p in R-'^'*"-'^, one may assume that, for all {x, t) & D 

p{x,t) = p{p) + Vp(p) ■ {{x tf-p} + 0{\{x tf-p\^) 

where p{p){^ 0) and Vp(p) are the corresponding values of the C^-extension of the original 
P- 

Then one has the expression 

exp |— Vc^ + Irp • a;(c)| exp |— icrp • ^ ^ ^1 ^ j j v{x, t)p{x, t)dxdt 



where 



= 7(t) + II{t) + III{t) 
I{t) 



and 



P{P) 

cxp |Vc2 + lr{ 1^ ^ j - p} • a;(c)| exp |icr{ ^ ^ j - p} • 

II{t) 

exp jVc^ + 1t{ ^ ^ j - p} • a;(c)| exp |icT{ ^ ^ j - p} 

xVp(p) • {(a; t)'^ — p}dxdt, 
III{t) 



-2cr ) f - ^ ^ 



-2cr ) V ^ 



exp|\/c2 + lr{ || j - p} • (^(c)| exp | icr{ f ) - p} • f 1 V 1 - 



t y \^ -2cT y V cV 

T „|2\ 



xO{\{xty -pY)dxdt. 

Here we prove that, as r — > oo I(t) ~ Cr""' with C 7^ and //(r) = III{r) — 0(t~^). 
First we study the asymptotic behaviour of /(r) as r — >■ 00. 
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Set 



Since 



/;ee-V^exp{.^e(.-p)-(_L)/^}cie 



r^5/(r)/p(p) = / . ^ ^^^^^ + 0(e-^^). (4.11) 

Using a parameterization of the segment Q, we have 
Lemma 4.4. T/ie formulae 

r dSjy) ^ ^ \yi - yo\ 

Jq (^s/^TT - iB{y, r))' {V^+l - iB{y,,T)) {V^+l - iB{yo, r)) 

and 

I 2Vpip).{y-p) ^ _ 

^( Vp(p)-(yi-p) ^_ Vp(p)-(yo-P) ^ 
(V^Mn:-iS(yi,T)) (V^Mn:-iS(yo,T)) 



we obtain 



(V?Tl - r)) (V?TT - iB{y^, t)) ' 



are valid. 



Proof. Let y = y{r]) be the parameterization of Q with y(0) = yo, = Z/i and 

y'(?7) = yi - yo- Since 

|^S(y,T) = S(yi,T)-S(yo,T), 

we have 

d ( 1 ^ _ . S(yi,r)-S(y2,T) 



dr] \y/c'' + l-iB(y,T)J (^^2TT - iS(y, r)) 

This gives (4.12) and a combination of a similar identity and integration by parts yields 
also (4.13). 

□ 

Since 

2c%i-p| r r 



S(yi,T) = ^^^rWl 
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and 



Biyo,r) = 



from (4.12) we have, as r — > oo 

|2 



/ 

Jc 



dS(y) . 5 \yi -yo\ 

' —1- 



^ - iB{y, r)) 2c _ ( ^ W 



A combination of this and (4.11) gives the formula: 

hm 2(cr)V/(r) = /'^^ " ^"''^^f^^ ^(^ 0). 



-p| ^Vc2 + 1 + 



\yi -P\ I Vc^ + l + i—\yQ-p\ j 

For the estimation of II (r) one has to make use of the growing factor in the 
imaginary part of the phase function in the integrand. For the purpose write 

d^T'^ JQ Jo 

Since 

■^0 (y^Tl-iB{y,T)) 



we have 



2Vp(p) • {y-p)dS{y) ^ ^^^^^^2g-r5V?+T) 



^ (Vc2 + l-iS(y,T)) 

Here, from (4.13) we have, as r — > oo 

2Vp{p) ■ {y-p)dS{y) iS \yi - yo\ Vp{p) ■ {yo - p) 



Q {V^Tl - tB{y, r)) 2c^ \y, - p\ /^___ ^ 
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iV 



This gives the estimate //(r) = 0{t ^) as r — > oo. The estimation of III{t) is rather 
easier than that of II{t): 

\in{r)\ < C f^ds jj?^)e-^rVc^\?^f\y-p\^dS{y) 
< rdre-^^'di I \y-p?dS{y) = 0{t-% 

Jo JQ 

Summing up, we conclude the existence of the nonzero limit of the integral 
T^exp |— -\/c^ + Irp ■ co{c)^ exp |— icrp ■ ^ ^ ^1 j-| J v{x,t)p{x,t)dxdt 
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as T 



oo. This completes the proof of Theorem 2.3. 
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